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We show that the nontwist phenomena previously observed in Hamiltonian systems exist also 
in time-reversible non-Hamiltonian systems. In particular, we study the two standard colli- 
sion/reconnection scenarios and we compute the parameter space breakup diagram of the shearless 
torus. Besides the Hamiltonian routes, the breakup may occur due to the onset of attractors. We 
study these phenomena in coupled phase oscillators and in non-area-preserving maps. 
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I. INTRODUCTION 

Dynamical systems are usually divided in two classes, 
conservative and dissipative systems. Although this divi- 
sion is conceptually useful, it can be misleading as far as 
there are systems that display both dissipative and con- 
servative (quasi-Hamiltonian) dynamics, not only for dif- 
ferent control parameters but also coexisting in different 
regions of the phase space. Systems with this counter- 
intuitive property are time-reversible [J but not Hamil- 
tonian. Two decades ago Politi et al. reported on a 
system of this type in a set of differential equations mod- 
eling a laser, and later other examples in arrays of Joseph- 
son junctions 0,3 and coupled phase oscillators |f| have 
been investigated. 

Recently, many fundamental studies in Hamiltonian 
dynamics have focused on nontwist systems, i.e., systems 
where the twist condition is violated. The twist con- 
dition asserts the nondegeneracy of the frequencies in 
the integrable regime and it is assumed in many fun- 
damental mathematical results [||. In generic situations 
we expect it to be locally violated. In such cases it was 
shown that genuine nontwist phenomena appear 0, • 
More recently, it has been shown that such nontwist 
effects are fundamental for the understanding of many 
different physical systems (see Ref. and references 
therein), e.g. , th e magnetic field lines in reversed shear 
tokamaks [jj, ll(& fill and the zonal flows in geophysical 
fluid dynamics |12| . 

In this paper we provide a novel link between these 
two aspects of classical dynamics by showing that the 
nontwist phenomena, studied so far in Hamiltonian 
systems, occur also in time-reversible non-Hamiltonian 
systems. In particular, we show the standard colli- 
sion/reconnections around the shearless torus and we 
study its breakup in a two-dimensional parameter space. 
Regarding the breakup of the shearless torus, we show 
that the usual Hamiltonian routes as well as a "dissipa- 
tive route" exist in time-reversible non-Hamiltonian sys- 
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terns. 

The paper is organized as follows. In Scc.[H]we present 
the definition of the twist condition and discuss the class 
of systems where it is expected to fail. In Secs. lIIII andllVI 
we present examples of the nontwist behavior in time- 
reversible non-Hamiltonian systems in continuous and 
discrete time, respectively. We summarize our conclu- 
sions in Sec. [V] 



II. TWIST CONDITION IN TIME-REVERSIBLE 
SYSTEMS 

The study of near-integrable Hamiltonian systems has 
led to one of the greatest successes of modern classical 
mechanics. Near-integrable Hamiltonian systems can be 
written as #(1,0) = H (I) + effi (1, 0) , where (1,0) 
are the action-angle variables of the integrable Hamil- 
tonian Hq. In this approach, many important mathe- 
matical results (e.g., the KAM theorem) are valid as- 
suming some nondegeneracy condition of the frequencies 
9 k = dH (I)/dI k (see, e.g., Appendix 8 of Ref. @). The 
most common nondegeneracy condition is the twist con- 
dition, defined as 
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for continuous and discrete time systems respectively. 
While the continuous version needs some modification 
when the dynamics is restricted to a specific energy shell 
||J, the discrete version can be directly applied to both 
maps and the reduced dynamics obtained by a Poincare 
section of the flow. In the (near-)integrable regime the 
twist condition assures that the frequency of the invari- 
ant tori varies monotonically with the action. Conversely, 
a local violation of the twist condition usually implies 
the existence of a torus with maximum or minimum fre- 
quency, the so-called shearless torus. Around the shear- 
less torus different nontwist phenomena were discovered 
in Hamiltonian systems, e.g., separatrix reconnections 
and island chains collisions [3, manifold reconnections 
of hyperbolic points in the chaotic regime |13| . mean- 
dering 0| and the fractality of the shearless torus at 
criticality @, 0, . These phenomena were observed 
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m area preserving maps and Hamiltonian flows 

with one and a half \VA and two [n| degrees of free- 
dom and they typically have a strong impact in transport 
properties of the system. 

In this paper we show the existence of the nontwist 
phenomena in time-reversible non-Hamiltonian maps and 
flows. A dynamical system is called time-reversible if 
there is an involution G (i.e., G 2 = Id) that reverses the 
direction of time. For example for dynamical systems de- 
scribed by a first order differential equation dx/dt = F(x) 
or by a mapping x n+1 = Lx„ reversibility implies 

^l = -F(G(x)) and LoGx„ +1 =Gx„, (2) 
at 

respectively. This condition alone does not ensure the 
observation of quasi-Hamiltonian dynamics. Addition- 
ally, the dimension of the invariant set of G should 
be large enough compared to the dimension D of the 
phase space. As discussed by Topaj and Pikovsky Q 
dim(Fix(G)) >D/2{D even), (D - l)/2 (D odd) is usu- 
ally required for continuous systems. 

The twist condition is defined for time-reversible 
non-Hamiltonian systems locally in the regions of the 
phase space where quasi-Hamiltonian dynamics is ob- 
served. Similar to the Hamiltonian case, a regime 
where the dynamics is locally integrable is needed in or- 
der to define action-angle variables (1, 9) to be used in 
Eq. (JTJ. Many fundamental results have been obtained 
for the quasi-Hamiltonian dynamics of time-reversible 
non-Hamiltonian systems Q (including the extension of 
the KAM theorem with the same optimal nondegeneracy 
conditions [20|). but until now there was no systematic 
study of the breakdown of the twist condition in these 
systems. 

A major question in this context is: in which 
kind of systems (Hamiltonian or time-reversible non- 
Hamiltonian) the violation of the twist condition is to 
be expected? Though we do not intend to give a general 
answer to this question, we note that systems possess- 
ing a symmetry in the phase space, apart from time re- 
versibility, exhibit naturally a shearless torus. Consider- 
ing that two tori related by this symmetry have the same 
frequency, and assuming continuity, we conclude that in 
the integrable limit there is at least one shearless torus 
between the two symmetric tori. This torus is invariant 
under the symmetry, what can be used (together with the 
involution G) to locate indicator points (IPs) [T(| . IPs are 
points that belong to the shearless torus whenever it is 
not broken. Further examples are dynamical systems de- 
scribed by phase variables (i.e., evolving on T D ), where 
the symmetry is given by the periodicity <j> — + 2n 
of the variables. This case is exemplified in Fig. 2] for 
a non-Hamiltonian time-reversible system composed of 
phase oscillators, described by Eq. JSJ below. Fig. QJa) 
shows that the phase space of the system is foliated by 
tori. Fig. ^b) shows the rotation number of these tori as 
a function of one phase space variable for three different 
control parameters s. Due to the periodicity of the phase 



variables, at least two shearless torus exist (maxima and 
minima of the frequency). As emphasized in Fig.^ one 
is located at the diagonal and the second contains the 
IPs tpi = -ijj 3 = ±7r/2. 
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FIG. 1: (Color online) (a) Poincare section (tp2 = tt/2) of 
the system in Eq. for ui = l,e = 0.2 where two nontwist 
tori are emphasized. IPs are marked with the symbol •. (b) 
Rotation number of the tori as a function of the coordinate 
ip3 at fixed tpi = — 7r/2 for ui = 1 and different values of e (see 
legend) . 



III. CONTINUOUS-TIME SYSTEMS 

Time-reversible non-Hamiltonian flows are encoun- 
tered in several physical situations (see for a survey) 
such as an externally injected class-B laser or arrays 
of N coupled differential equations |3|, |5j . One system 
where the properties discussed in Sec. [H] can be found 
is an array of N coupled phase oscillators like the one 
considered by Topaj and Pikovsky [j| 

<Pk = ^k+sfifk-i - fk) + ef(<Pk+i - <Pk), k=l, ...,N 

(3) 

with boundary conditions ifQ = tpi, <Pn+i = VJV and 
vanishing coupling when the phases of the oscillators are 
equal, i.e., /(0) = 0. Taking the phase differences ip k = 
fk+i — fk the number of variables is reduced by one: 

4k = A fc +e/(V>fc_i)+e/(^fc+i)-2e/(V'fc), * = 1, N-l 

(4) 

with A k = fl k+ i - ft k . 

If / is an odd function and the natural frequencies 
are taken symmetrically = Ajv-fe the system is re- 
versible. The associated involution is G : ipk — ► 
7r - Tjj N - k ; and Fix(G) is given by ipk + 4>N-k = tt. In 
order to visualize the nontwist phenomena mentioned in 
Sec. II we restrict ourselves to three variables [N = 4 
in Eq. , and we analyze the dynamics by means of a 
two-dimensional Poincare section at ip2 — tt/2. The set 
of differential equations [for concreteness we already take 
/(•) = sin(-)] we will study is: 

ipi = uj — 2esin^i + esin^, 

ip2 = 1 — 2esin^2 + ssin^i + esin^3, (5) 

-03 = u — 2e sin 1P3+ s sin tp2 , 
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where (rescaling coupling and time) we may set A2 = 1. 
In addition to the coupling variable e, we take u) = Ai = 
A3 as a second parameter of our system. This allows a 
better exploration of the behaviors of the systems (e.g., 
codimcnsion-two points); and it is also useful in order 
to compare our results to those obtained in the so-called 
standard nontwist map by two-parameter sweep [bH Il6| . 

Figure^a) shows a Poincare section of the phase space 
for small coupling e. Phase space is foliated by tori with 
different rotation numbers r [Fig. IHb)]. Equation (JSJ is 
invariant under the transformation Tpi 3 — > i which 
imposes the existence of the invariant manifold ipi — ip3 
that corresponds to the diagonal shearless torus in Fig.Q] 
The symmetry with respect to this diagonal torus to- 
gether with the involution G that maps a torus into 
itself allow us to locate the IPs (ipi = — V>3 = ±tt/2) 
for the off-diagonal shearless torus showed in Fig. [we 
will refer to it as the shearless torus (ST)]. For the case 
e = 0.2 illustrated in Fig.^a) there are no further shear- 
less tori. However, increasing the control parameter the 
rotation number of these tori may pass from maximum 
to minimum (or vice- versa) giving birth to other pairs of 
shearless tori [e.g., e = 0.3 shown in Fig. Gfb)]. We are 
specially interested in the region around the ST, where 
we will report in what follows the nontwist phenomena 
mentioned in Sec. II. 



A. Collision/Reconnection Scenarios 

We start describing the two standard colli- 
sion/rcconncction scenarios studied so far in nontwist 
Hamiltonian systems Q • Figure [2] shows the two sce- 
narios for our system: periodic orbit collision (a-c) and 
separatrix reconnection (d-f). We find the same phe- 
nomena as in Hamiltonian systems, the only appreciable 
difference being the local non- volume preservation of the 
flow J^}. Since in the quasi-Hamiltonian dynamics the 
phase-space volume is preserved only in time-average, 
we see in Figs, ^a) and [21 bunches of tori compressed 
and expanded in different regions of the phase space and 
Poincare-Birkoff chains composed of islands of different 
sizes. Apart from these nonsymplectic features, Fig. [21 
essentially displays the two standard scenarios: 

(i) Figure [2fa-c) shows the collision sequence of two 
symmetry related Poincare-Birkhoff chains with 3:4 
rotation number, by increasing lu at fixed e. The 
collision occurs in Fig. 2(b), and in Fig. 2(c) one 
may see the resulting dipolar structures formed by 
two saddle cycles on the torus and one center at 
each side of the ST. When ui is increased further 
the dipolar structure shrinks and finally disappears 
(at the point where the two saddles cycles on the 
torus annihilate each other). 

(ii) Figure GId-e) shows the reconnection scenario. 
Twin Poincare-Birkhoff chains with 2:3 rotation 
number are out of phase [Fig. [2{d)] and approach 



each other with the ST between them. At a critical 
parameter [28| the separatrices of both islands re- 
connect [Fig.[2Je)]. Finally, passed the bifurcation, 
the torus exhibits a characteristic meandering, due 
to the exchange of saddles between the two involved 
chains of islands. 



Scenario (i) 
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Scenario (ii) 




FIG. 2: (Color online) Poincare section of system (J^J for fixed 
e = 0.25 and different values of lu. Sequence (a) lu = 0.868, (b) 
lu — 0.8687606, and (c) lu — 0.869 illustrates the collision of 
3:4 island chains. Sequence (d) lu = 0.801, (e) lu = 0.801523, 
and (f) lu = 0.802 illustrate a reconnection around 2:3 reso- 
nances. 

As usually found in the literature, scenario (i) is ob- 
served for rotation numbers with even denominators, 
whereas scenario (ii) applies to odd denominators. Nev- 
ertheless, as demonstrated in [ll|, scenario (i) is not ex- 
pected in the absence of symmetries [scenario (ii) is then 
observed for both even and odd cases] (see, e.g., [ijj). 
As final remarks we mention that the shearless torus on 
the diagonal also exhibits the collision scenario [type (i)] . 
Heuristically the reconnection scenario [type (ii)] cannot 
occur because this torus is constrained to the diagonal 
and cannot meander. Another remark is that for reso- 
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nances with small denominator (e.g., 1:1 and 1:2) a more 
complicated reconnection scenario occur due to the in- 
terplay of simultaneous resonances of both shearless tori 
(the diagonal one and the ST). 



B. Destruction of the shearless torus 

Our next numerical experiment was to vary the two 
control parameters (e,uj), in order to find the region of 
parameter space where the ST exists, and in this way un- 
derstand the routes for the destruction of the ST. Due to 
the rather difficult numerical implementation of our sys- 
tem (if compared to a two dimensional map), we needed 
to find an efficient method to detect the existence of the 
ST for a large set of parameter values. By far, the most 
efficient way we found was to resort to the Slater's theo- 
rem [2l]]. In short, we start an initial condition at one IP 
and we check whether the number of different recurrence 
times to a small box around this point is at most three 
(see Appendix for details). 
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region between the two 1:2 lines in Fig.|3J). As expected, 
the torus breakup is favored by the global bifurcations 
at the collisions/reconnections, and accordingly there is 
a match between the solid lines and the white "wedges" 
in Fig. |21 Singular points at the tips of the shaded region 
correspond to noble irrational rotation numbers where 
the torus is expected to be fractal 0,113, Fig.0{a,b) 
show the torus just before and after destruction. 
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FIG. 4: (Color online) Two different routes for the breakup of 
the shearless torus: (a)-(b) Hamiltonian-like through critical 
point and (c)-(d) dissipative. (a) torus near criticality e = 
0.425256, uj = 1.0335 (the inset shows the fractal structure 
of the torus), (b) Torus after criticality e = 0.425257, lu = 
1.0335 (the inset show that the torus is destroyed), (c) Near- 
integrable phase space e = 0.1, ui = 0.2. (d) Attracting fixed 
point e = 0.11,0/ = 0.2. 



FIG. 3: (Color online) Parameter space (e,ui) of Eq. 0. 
The shaded region indicates the parameters where the ST 
exists; the inset evidences the fractal structure of the border. 
Several solid lines indicate the collision/reconnection param- 
eters for rotation numbers 1:2, 2:3, 3:4 and 1:1. The symbols 
■ a, A, ♦) are located at parameter values with characteris- 
tic dynamics shown in the corresponding figures. The dotted 
line FP marks the onset of an attracting fixed point in the 
system. 

The result of our calculation is shown in Fig. [3] The 
shaded region indicates the parameter values where the 
ST exists, while in the white zone the ST does not exist. 
The fractal character of the border limiting these two re- 
gions is apparent (see inset) (cf. 0,0). The solid lines 
indicate the regions where a few collisions (1:2, 3:4) and 
reconnections (2:3,1:1) occur. Regions where a couple of 
saddle orbits lie on the ST [as in Fig. Etc)] are consid- 
ered as regions where the ST exists (hence the shaded 



Since the system is not Hamiltonian the existence of 
attractors cannot be excluded |5(. For instance, below 
the dotted line (denoted FP-line) in Fig. (0 a globally 
attracting fixed point exists (a periodic orbit in the full 
phase space) . Dissipation introduces alternative ways for 
the destruction of the ST, as shown in Fig. 0tc,d). The 
bifurcation mechanism for the system J5j follows: At the 
FP-line a marginally stable fixed point is born on the 
invariant set of G at tpi = 7r/2 (= ^3). This point be- 
longs to the diagonal shearless torus, which is broken 
only through this bifurcation. Below the FP-line four 
fixed points exist: one attractor, one repeller (mirror of 
the attractor as demanded by reversibility) , and two sad- 
dles on the invariant set of G [see Fig.^Jd)]. Going back 
to the original system of N = 4 phase oscillators, FP- 
line marks the onset of synchronization of oscillators 1-2, 
and 3-4 (i.e. a two-cluster state). Note that the usual 
saddle-node bifurcation in this kind of transition to syn- 
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chronization is forbidden due to reversibility: for every 
attractor a mirror repeller must exist. In a wide region 
above of the FP-line a (chaotic) attractor and a set of 
Hamiltonian-like tori surrounding the diagonal torus co- 
exist. 



C. Other systems 

In addition to the system studied above, Eq. (J5J, it is 
interesting to mention the case of an array of N Joseph- 
son junctions subject to a parallel resistive load. This 
case may be described by globally coupled phase oscilla- 
tors (see Q where the case N — 2 is studied in detail). 
For N = 3 and assuming that two of the junctions (say 
j = 1 and j = 3) are identical, we get a set of equations 
with the same symmetries of Eq. ||SJ| 



6i = ft + a sin (j>i 
t>2 = f2 + a sin 4>2 
63 — £1 + a sin </>3 



3-Ej=l sin( ^ 
1 \p3 
3 2.^1 



(6) 



We have observed in this system the same colli- 
sion / reconnection scenarios described previously for sys- 
tem J5J. As before, the diagonal torus 4>\ = ^3 ceases to 
exist by the onset of a stable fixed point (of the Poincare 
section tp2 = tt/2). However, differently from system (JSJ, 
at first this stable fixed point is not globally attracting 
(as in the N — 2 case studied in 3]). Thus, the "dissi- 
pative breakup" of the off-diagonal shearless torus does 
not occur right after the bifurcation but only when the 
basin of attraction of the fixed point becomes the whole 
phase space. 



IV. DISCRETE-TIME SYSTEMS 

A simple procedure to obtain a time-reversible map is 
to integrate a time-reversible differential equation x = 
F(x) semi- implicitly |29j : 



x n+ i - x„ = fc[F(x„ + i) - F(x„)], 



(7) 



where k measures the finite-time integration step. Ap- 
plying this procedure to equations of phase oscillators, 
similar to those studied in Sec. 11111 we have obtained 
time-reversal maps on the torus where the same non- 
twist phenomena were observed. Certainly, the fact that 
mapping Q is defined implicitly has several drawbacks, 
however, avoiding details, we note that the methods to 
construct (explicit) time-reversal mappings described in 
Ref. are not suited for the torus topology. In this sec- 
tion we study a different two-dimensional time-reversible 
explicit map defined on a cylinder, which illustrates that 
nontwist phenomena are not restricted to dynamics on a 
torus. 



Consider the following map L 

y n + asin(27ra n ) 



Vn+l 
X n +1 



1 + by n sva(2'Kx n 
x n + cos(27ry n+ i) 



mod(l) 



(8) 



which is a particular example of a class of time-reversible 
maps discussed on page 103 of Ref. pj. Equation © 
diverges for 1 4- by n sin(27r:E n ) = and thus we restrict it 
to the invariant region — y? < V < (limited by the 
invariant lines y = ±^/g) and the control parameters to 
0<ab<l. 

The map (JSJ) can be written as L = Mi o M% , where 
Mi, 2 are the following involutions 



Mi 
M 2 



x' = -x,y' 



y + asin(27ra;) 

1 + bysin(27rx) ' 
cos(27ry),y' = y. 



This property ensures time reversibility under G 



M 2 
J = 



1].' The determinant of the Jacobian is given by 
i-absm (271-x) ^ shows that the map is Hamilto- 
nian only in the case 1) = (when one recovers the Harper 
map [Xjj)- In the limit of small control parameters the 
dynamics of map JSJ) can be considered as a perturba- 
tion of an integrable Hamiltonian system and it remains 
integrable almost everywhere since the KAM theorem, 
generalized to time-reversible systems, applies pflj . 




FIG. 5: (Color online) Parameter space of the map |0J: 
shaded regions indicate the parameters where the ST ex- 
ists. Above the curve b = 1/a the map is not defined. 
Lower inset: magnification showing the white wedge struc- 
ture in finer scale. Upper inset: phase space for parame- 
ters a = 0.11, b = 5. 



In Sec. IIIII we argued that when at least one symmetry 
is present in the system, apart from time reversibility, the 
twist condition must be violated in the near-integrable 
regime. In the case of map (jHJ this symmetry is given 
by x 1 — x + i, y' = —y, which can be used (together with 
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M 2 ) to determine the indicator points x — 0.25, y = 
and x — 0.75,?/ = 0. According to Eq. the twist 
condition is violated in map |H| at y = — asin(27ra;). In 
the upper inset of Fig. [3] we show the phase space of the 
map © for a = 0.11, b = 5. The ST intersects the 
curve where the twist condition is violated Q . Using the 
method described in the Appendix we determined the 
regions of the parameter space (a, b) where this torus ex- 
ists. The results are shown in Fig. [5] where a fractal- like 
border similar to the one observed in Fig. [3] is clearly 
recognizable. Indeed, we have observed numerically that 
all nontwist phenomena discussed in Sec. IIIII arc also 
observed in the time-reversible non-area-preserving map 
given by Eq. JHJ. 



deed conservative dynamics has been already observed 
in phase models such as the dynamics of solitary waves 
in chains of dispersively coupled oscillators 0, and the 
finite-dimensional Kuramoto model (25|. The role of 
time-reversibility and the violation of the twist condition 
in these non-Hamiltonian models remains to be under- 
stood. 
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V. CONCLUSIONS 

Timc-reversibility does not appear exclusively in 
Hamiltonian systems. Remarkably, there are time- 
reversible non-Hamiltonian systems that may exhibit, in 
addition to dynamics approaching an attractor, quasi- 
Hamiltonian dynamics. While major results in this con- 
text have been achieved in the past twenty years, a com- 
plete understanding of which features of Hamiltonian dy- 
namics can be extended to reversible systems is still lack- 
ing. 

Frequently, time-reversal symmetry is not the only 
symmetry of a system, in which case we have argued that 
the twist condition (i.e., the nondegeneracy of the fre- 
quencies) is violated. We have investigated the existence 
of nontwist phenomena, previously studied in Hamilto- 
nian systems, in time-reversible non-Hamiltonian sys- 
tems. From our results, obtained in both continuous- and 
discrete-time systems, we may conclude that the same 
nontwist phenomena are reproduced in non-Hamiltonian 
systems. In order to study the breakup of the shear- 
less torus we have developed a novel numerical method 
to compute the parameter space breakup diagram (see 
Appendix). Our method is potentially useful in a wide 
class of problems concerning the detection of quasiperi- 
odic motion in a multiparameter space. We have iden- 
tified the parameter regions of the usual Hamiltonian 
routes of breakup of the shearless torus (through col- 
lision/reconnection and through a critical- fractal torus) 
and the parameters where the breakup is due to the onset 
of attractors. 

Besides the theoretical interest of expanding the class 
of systems where the nontwist phenomena occur, our re- 
sults are relevant for the comprehension of specific time- 
reversible non-Hamiltonian systems. Many dissipative 
systems posses an attractor of limit-cycle type where a 
marginally stable phase dynamics exists. The interaction 
of such systems may lead to quasi-Hamiltonian phase dy- 
namics, in which case nontwist phenomena should be ex- 
pected. Usually a large number of oscillators are coupled, 
what provides additional motivation for the extension of 
the nontwist phenomenology to higher dimensions. In- 



Appendix 

The knowledge of the critical parameter values of a 
map where an invariant rotational circle breaks can be 
crucial in many situations. We briefly describe here the 
techniques previously employed to calculate the breakup 
diagram of the shearless torus and we introduce a new 
method which has proven to be computationally efficient. 
Even though we have applied our method to compute the 
breakup diagram it is rather general and we believe that 
it can be useful whenever quasi-periodic motion has to 
be efficiently detected. 

The first breakup diagrams were calculated 0| con- 
sidering the torus to be destroyed whenever a trajectory 
started at an IP leaves a region that certainly would con- 
tain the curve. A slightly different method was proposed 
in Ref. [14( using the fact that a trajectory in a periodic 
or quasi-periodic motion leads to a converging winding 
number. These procedures can be computationally ex- 
pensive: in a typical mixed phase space we expect that a 
trajectory will stick around the complex structure of can- 
tori and regular islands while wandering into the chaotic 
sea and the time needed to detect the non-existence of 
the tori could be very large. A refined method relies on 
Greene's residue criterium and can be used when high ac- 
curacy is needed (e.g., to precisely locate critical points 
in the breakup diagram corresponding to noble rotation 
numbers) while it is not suited to explore a large param- 
eter region |sH22||. 

Here we take advantage of a simple property of ro- 
tations in one dimension to develop a method that is 
together fast and quite robust, and hence suited to scan 
a large set of parameters. The quasi-periodic dynamics 
restricted to the torus can be reduced to a simple rota- 
tion of the circle using a natural parameterization of the 
curve (an example can be found in |26j). Slater's theo- 
rem [21| states that for any irrational rotation uj and for 
any connected interval there are at most three different 
return times. Moreover, in the case of three different re- 
turn times one of them is the sum of the other two and 
two of them are always consecutive denominators in the 
continued fraction expansion of the irrational number lu. 
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Accordingly, our method simply consists in counting 
the number of different return times of the iterates of the 
IP xq inside an arbitrary region that contains a connected 
part of the torus around xo (for the specific cases Figs. [31 
and [5] we used a box centered at one IP [3Cj). The torus 
is considered to be broken whenever the different return 
times violate the conditions imposed by Slater's theorem 
(e.g., whenever their number exceeds three). When the 



rotation number of the torus is not known a priori (as in 
our case) the additional restriction that uses its continued 
fraction expansion can not be used. 

The implementation of the procedure is straightfor- 
ward provided one point on the torus is known (IP). It is 
evident that the constraints on the return times are quite 
restrictive and they are typically rapidly failing when the 
torus is broken and the trajectory enters the chaotic sea. 
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